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Abstract 

The timc-dcpcndcnt transmission coefficient for the Kramers problem exhibits dif- 
ferent behaviors in different parameter regimes. In the high friction regime it decays 
monotonically ("non-adiabatic"), and in the low friction regime it decays in an oscil- 
latory fashion ("energy-diffusion-limited"). The generalized Kramers problem with an 
exponential memory friction exhibits an additional oscillatory behavior in the high fric- 
tion regime ("caging"). In this paper we consider an oscillatory memory kernel, which 
can be associated with a model in which the reaction coordinate is linearly coupled 
to a nonreactive coordinate, which is in turn coupled to a heat bath. We recover the 
non-adiabatic and energy-diffusion-limited behaviors of the transmission coefficient in 
appropriate parameter regimes, and find that caging is not observed with an oscilla- 
tory memory kernel. Most interestingly, we identify a new regime in which the time- 
dependent transmission coefficient decays via a series of rather sharp steps followed 
by plateaus ( "stair- like" ) . We explain this regime and its dependence on the various 
parameters of the system. 



1 Introduction 



The classic Kramers formulation of reaction rates in solution M and its generalization 
to non-Markovian solvents 0] has provided many theoretical challenges over the past six 
decades [||-||])- In this formulation the reaction coordinate x(t) is modeled as evolving in a 
double-well potential V(x) with a barrier separating the reactant and product states. The 
solvent effects are modeled in terms of fluctuating and dissipative forces. A full understand- 
ing of the dependence of the rate coefficient k on the dissipation in the Markovian solvent 
limit (the "turnover problem") has only been achieved in the last few years |@,@, Com- 
parably thorough understanding in the case of a non-Markovian solvent is not yet available. 
Understanding of the temperature dependence of k is also far from complete |S-12|. Clearly, 
there is yet a great deal to learn about this classic problem. 

In the past two decades, and most especially in the past few years, attention has also 
been paid by a number of investigators to the time-dependence of the rate coefficient, that is, 
the way in which k(t) approaches its asymptotic value k(oo) [11-16]. This time dependence 



directly mirrors the dynamics of the reaction coordinate in the barrier region on the way 
toward capture by one well or the other. Our focus is on this time dependence and the way 
that it is influenced by the parameters of the system. 

The generalized Kramers problem is based on the dynamical equations for the reaction 
coordinate 

* dt 1 T(t - x(f) - dVef J {x) + F(t), (1) 



o 



dx 



where a dot denotes a time derivative, V e ff(x) is an effective potential related to V(x) (cf. 
next section), T(t — t') is the dissipative memory kernel (which we will often simply call 
the memory kernel), and F(t) represents Gaussian fluctuations that satisfy the fluctuation- 
dissipation relation 

(F(t)F(t')) = k B TT(t-t>). (2) 

The brackets (• • • ) denote an ensemble average, ks is Boltzmann's constant, and T is the 
temperature. The fluctuations and dissipation account for the interaction of the reaction 
coordinate with the surrounding medium. The original Kramers problem dealt with a 
Markovian solvent, that is, with instantaneous dissipation: 

r(t) = 2 7 a(t). (3) 

The parameter 7 is the dissipation parameter or damping parameter. Generalizations to 
the non-Markovian problem have typically focused on exponential memory kernels 0, 

T{t) = 1 -e- t '\ (4) 
r 

and on oscillatory memory kernels Q, 

r(t) = r(0) e- t/r ^cos at + ^- sin Sltj (5) 

(this memory kernel and the parameters in it will be discussed in detail in the next section) . 
Another generalization, which we do not address in our work, deals with Gaussian memory 
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kernels @,|T5|, 



r (i ) = (^ 1/2 V^. (6) 



7T / 

In all of these generalizations r is a measure of the decay time of the memory kernel or, 
equivalently, of the correlation time of the fluctuations. 

In subsequent sections we will provide a brief graphic review of the results addressed 
in our previous work, which succinctly are as follows. The time-dependent rate coefficient 
for the Markovian solvent at high damping (that is, beyond the "turnover" regime) decays 



monotonically towards its equilibrium value [16|. For the exponential memory at high 
damping there are two distinct types of time dependences, the "non-adiabatic" , in which 
the rate coefficient decays monotonically to its equilibrium value (as in the Markovian 
case), and the "caging", in which the decay to equilibrium is oscillatory with a frequency 



characteristic of an effective caging potential [16|. At very low damping (that is, below the 



"turnover"), the decay of the rate coefficient to its equilibrium value is again oscillatory, 
but now with a frequency characteristic of the bistable potential. This behavior is apparent 
for the Markovian solvent [11| and also for the exponential memory p) in this energy- 



diffusion-limited regime. We have shown that the theoretical predictions agree very well 
with numerical simulations for all of these generic behaviors. 

In this paper we complete our analysis with a study of the oscillatory memory kernel. 
In appropriate limits we recover the typical low-damping behavior of the rate coefficient 
and also the high-damping non-adiabatic monotonic behavior, although caging, as we will 
see, can not be achieved with an oscillatory memory. Most interesting, perhaps, is the ap- 
pearance of a new time dependence different from those previously observed or anticipated. 
This new time dependence is a consequence of the new features of the memory kernel such 
as the fact that it alternates between positive and negative values. We will present and 
explain this new behavior, and determine the parameter regimes where it may be observed. 

This paper is organized as follows. In Sec. Q we introduce the model Eq. (ffl) in detail. 
In fact, we present two equivalent versions of the model. One version invokes a solvent 
coordinate which is coupled to the reaction coordinate and also coupled to a heat bath. This 
double presentation not only clarifies the physical origin of the oscillatory memory kernel, 
but it also leads to more transparent interpretation of the resulting time dependence of the 
rate coefficient. In Sec. ||| we describe our simulations and numerical procedures. Section ||| 
presents a graphical summary of the various time dependences obtained numerically in 
earlier work and provides a context for the presentation of the new behavior identified in 
the oscillatory memory system. In Sec. || we discuss analytic approximations that serve as a 
backdrop for our analysis and detailed explanation of the new behavior, which is presented 
in Sec. ^. The results and conclusions are summarized in Sec. [/]. 

2 The Model 

We first present an alternative (two-variable) model that eventually leads to Eq. (p]) with 
Eq. (H) . The potential energy of the two- variable model is 

V(x,y) = V(x) + ^-y 2 + ^(x-y) 2 , (7) 
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where the solvent is explicitly represented by a harmonic coordinate y. The reaction coor- 
dinate x evolves in a bistable potential that is taken to be of the familiar form 



V{x) = ^-lf. (8) 

We take Vq as the energy unit throughout this work and thus set it equal to unity. The 
reaction coordinate is coupled to the solvent coordinate via a harmonic spring of force 
constant k. The dynamical equations for the coupled system, assuming that y is coupled 
to a heat bath at temperature T, are 

dV(x) 

x = h k(y - x), 

ax 

V = ~u 2 y + k(x -y) -jy + f(t). (9) 

Here 7 is the friction coefficient for the solvent coordinate and f(t) represents 5-correlated 
Gaussian fluctuations that satisfy the fluctuation-dissipation relation 

(f(t)f(t')) =2 1 k B T8{t-t'). (10) 

Throughout we take the barrier height to be large compared to the temperature, ksT < 1/4. 
We call Eq. (g) the "extended" representation of our system. 

Although the extended model can readily be integrated numerically, initial conditions 
are not arbitrary and require careful consideration. The reduction of the model (JoJ) to 
the generalized Kramers problem (||) with the fluctuation-dissipation relation @ requires 
distributions for the initial solvent coordinate y(0) and velocity y(0) that satisfy certain 
conditions (cf. Appendix |A|, where these initial conditions are presented in detail). 

A number of theoretical approaches to this problem deal, instead, with the completely 
equivalent "contracted" or "reduced" representation obtained by explicitly integrating out 
the solvent coordinate y. Among these is the work of Grote and Hynes || and that of 
Kohen and Tannor (KT) [p!q| . The resulting equivalent single- variable problem is shown in 
Appendix^] to be given by Eq. (ffl) with the effective potential 

V eff {x)=V{x) + \-^ l x\ (11) 

Depending on the relative values of parameters, the resulting memory kernel can decay 
monotonically ( "hyperbolic" case) or it can decay in an oscillatory fashion ( "trigonometric" 
case). We are specifically interested in the trigonometric case, which requires that 

-J 2 - k <0. (12) 
The associated memory kernel is 

Tit) = e~ 2* (-3- sin Sit + cos Sit) (13) 

u> + k \2\l ) 
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with the frequency 



^yWfc-g). (14) 

We expect that the oscillatory character of the memory kernel may lead to new regimes of 
dynamical behavior that will become evident in the time dependence of the rate coefficient. 
We do not pursue the hyperbolic case because we expect behavior similar to that found 
earlier for the exponential memory kernel and hence do not expect new behaviors in this 
case. 
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Figure 1: Solid line: original bistable potential V(x) of Eq. (||). Dotted line: effective 
potential V e ff(x) of Eq. ( |i~l| ) for J 2 = 0.5 and k = 1. 



From the explicit form (|ll| ) of the effective potential we see that the additional quadratic 
term moves the minima of the wells of the bistable potential V{x) from ±1 to 



Ml 



0J 2 k 

J 2 + k 



(15) 



and diminishes the barrier from 1/4 to 



AV 



1 1 



eff 



UJ 2 k 



2 \u 2 + k 



+ 



1 ( UJ l k 

4 V + k 



(16) 



Both effects can be seen in Fig. ffl. It is easily shown that the barrier disappears entirely 
when uj 2 k > lo 2 + k, at which point the very nature of the problem changes. We thus 
constrain our parameters k and uj 2 to ensure bistability: 



u 2 k < u 2 + k. 



(17) 



To summarize, then, the model to be considered in this paper is given by Eqs. (EI) with 
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the potential (||) and the fluctuation-dissipation relation (|l0|) (extended representation), or, 
completely equivalently, by Eq. (||) with the effective potential (^TJ), the memory kernel 
di~3|), and the fluctuation-dissipation relation (|2|) (reduced representation). Whichever for- 
mulation is used, our parameters are constrained by the inequality (0), which ensures a 
bistable effective potential, and by the inequality (|l2|), which ensures an oscillatory memory 
kernel. For some purposes the extended representation provides a more convenient view- 
point, while for others the reduced representation is more transparent. In particular, we 
find the extended representation more convenient for numerical simulations. 

As a final point in this section it is important to note the altered significance of param- 
eters in the oscillatory memory kernel compared to the Markovian or exponential models. 
In the latter two cases "high friction" and "low friction" refer to the value of 7 since this 
parameter directly measures the strength of the dissipative force that extracts energy from 
the reaction coordinate into the bath. Thus, in the exponential memory case 7 is the 
value of r(0) and also of the integrated memory kernel. In the case under consideration 
here, however, 7 is a measure of the dissipative force on the solvent coordinate and not 
directly on the reaction coordinate. Although 7 indirectly affects the loss of energy of the 
reaction coordinate, the energy loss channel is now principally determined by the coupling 
strength between the reaction coordinate and the solvent coordinate. Correspondingly, now 
r(0) = k 2 /(uj 2 + k). The integrated memory kernel is r(0)7/(u; 2 + k), thus reflecting the 
overall influence of 7. However, it is the coupling constant k that now essentially deter- 
mines whether we are in the "high friction" or "low friction" regime (more details will be 
presented in Sec. |j). 



3 Simulation Method: Initial Conditions and Other Details 

The quantity of interest is the time-dependent rate coefficient k(t) for an ensemble of parti- 
cles evolving according to x(t) in Eq. (|l|) or Eq. (^). Numerically we find it more convenient 
to work with the extended system The coefficient k(t) is the time-dependent mean rate 
of passage of the particles across the barrier at x = 0. The usual focus on the deviation of 
k(t) from its equilibrium transition state theory (TST) value leads to the expression 

k(t) = K{t)k TST (18) 

where k TST = (\^2/tt) exp(— l/AksT) is the transition state theory rate that assumes that 
particles never recross the barrier. The transmission coefficient k{€) is the correction to 
transition state theory that includes both the temporal dynamics and the effects of those 
particles that do recross the barrier. We are interested in the dynamics of the transmission 
coefficient K,{t). 

Numerically, one might try to calculate kit) directly by starting all the particles in one 
well and computing at each time how many of them have crossed to the other well. It 
would require an exceedingly long calculation to gather statistically significant data in this 
manner, since the reaction barrier is very much higher than the thermal fluctuations. The 
reactive flux formalism |3,17, ij| that relies on Eq. (|i~8| ) overcomes this difficulty since the 



transmission coefficient can be calculated by dealing only with an ensemble of particles 
whose initial position is above the barrier [x(0) = 0]. The slow process of "getting there" is 
already included in k TST . Half of the particles that start above the barrier have a positive 
velocity distributed according to the Boltzmann distribution in energy, and the other half 
have the same distribution but with negative velocities. 
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Upon imposing the initial conditions on y discussed in Appendix [A|, we run quite a few 
iterations for the solvent coordinate evolution in order to obtain even better thermalization. 
Having achieved this, we then integrate the fully coupled system (H) with the following 
distributions for the initial reaction coordinate position x(0) = x Q and initial velocity x(0) = 

P(x )=5(x ), (19) 



PM-fe^f-jjjfrl (20) 



The numerical integration is carried out using the second order Heun's algorithm [IS, 2(J. 
In all our runs our ensemble consists of N = 10, 000 particles and we use a very small time 
step (At = 0.001). The transmission coefficient is calculated from these simulated data 
according to the relation [17] 

n(t) - Z±® - (21) 
K[t) -N + (0) iV_(0)' {21) 

where N+(t) and iV_(£) are the particles that started with positive velocities and negative 
velocities respectively and at time t are in or over the right hand well (i.e., the particles 
for which x(t) > 0). Alternatively and completely equivalently (via a simple symmetry 
argument) one can start all the x particles with a positive velocity and then 

K(t) =K+(t)-K-(t) (22) 

where n+(t) is the fraction of particles that are in or over the right hand well at time t, and 
ft_ (t) is the fraction in or over the left hand well. Furthermore, it is easily argued that |T| 

K_(t) = l-K+(t), (23) 

so it is sufficient to follow one or the other. 



4 Numerical Results 

The number of independent parameters in the generalized Kramers problem with an oscil- 
latory memory kernel is of course larger than for Markovian or exponential frictions. The 
values of n(t) and n st now in general depend on k, u 2 , 7, and ksT. Indeed, a systematic, 
even qualitative study of the transmission coefficient as was done, for example, in 1 1 1 . 0] 



would be quite complex. We focus on a more modest goal, that is, to capture qualitatively 
the different types of temporal behavior of n(t) and the broad parameter regimes where each 
occurs. These include the three regimes identified for the exponential memory, namely, the 
energy-diffusion-limited, the non-adiabatic, and the caging regimes, as well as possible new 
behaviors. 

The oscillatory memory kernel can exhibit different appearances depending on the pa- 
rameter choices. Figure [2] exhibits three distinct "generic" appearances, each roughly repre- 
sentative of a distinct parameter regime. Two of these mimic behaviors of the exponential 
memory kernel and might be expected to lead to transmission coefficients similar to those 
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Figure 2: Memory kernel F(t) vs t for the different regimes studied in Sec. |j. Solid curve: 
uj 2 = l.o, k = 0.14, and 7 = 0.667. Dashed curve: J 1 = 0.01, k = 0.75, and 7 = 1.74. 
Dotted curve: uj 2 = 0.01, k = 0.3, and 7 = 0.05. 



obtained earlier. The third, the strongly oscillatory kernel, is new and might be expected 
to lead to new behavior. Let us consider each case in turn, along with the resulting trans- 
mission coefficients. 

The solid curve kernel in Fig. ^ mimics the exponential memory kernel in the energy- 
diffusion-limited regime. The kernel T(t) is small at all times. In the case of the exponential 
memory kernel (Q) this behavior was insured by choosing 7 to be small, but in the oscillatory 
memory case, as noted earlier, the meaning of the parameters is different. Now a small 
value of r(i) and, in particular, a small value of T(0) is insured if we choose small values of 
k 2 /(u> 2 + k), that is, k must be small and/or uj 2 must be large. Note that the choices must 
still obey the constraint (|l~7|), but this is not a problem. The value of 7 is not constrained 



by the low friction requirement, but it is constrained by Eq. (12) if we want to insure that 
we are in the oscillatory regime. The value of 7 determines the oscillation frequency of the 
memory kernel T(t) but not its magnitude. 

Typical values of the parameters that satisfy the conditions to produce energy-diffusion- 
limited behavior while preserving the oscillatory character of the memory kernel are to 2 = 
1.0, k = 0.14, and 7 = 0.667. These are the values used to produce the solid curve in Fig. |2|. 
In this regime the low dissipation causes the dynamics of the system to be dominated by 
the slow variation of the energy and consequently by the repeated inertial recrossing of 
the barrier before the particles are trapped in one well or the other. As expected, we 
find the typical energy-diffusion-limited behavior for K,(t) (oscillatory curve in Fig. ||) that 
consists of a very small initial decay to a plateau up to a time beyond which K,(t) decays 
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Figure 3: Solid curves: the three typical behaviors of the transmission coefficient with 
oscillatory memory friction (numerical results). Oscillatory curve: ui 2 = 1.0, k = 0.14, 
and 7 = 0.667 (energy-diffusion-limited regime). Monotonic curve: uj 2 = 0.01, k = 0.75, 
and 7 = 1.74 (non-adiabatic regime). The temperature for these two cases is k^T = 0.025. 
Stepped curve: to 2 = 0.01, k = 0.3, and 7 = 0.05 (new "stair-like" regime). The temperature 
is ksT = 0.015. The dotted, dashed, and dot-dashed curves correspond to the KT theory 
results for the same parameter values. 



in an oscillatory manner to its equilibrium value. As shown and discussed in our earlier 
work 1 1% 12], the first decay is due to the few low-energy particles that immediately change 
their initial direction due to a thermal fluctuation and are trapped in the well opposite to 
the one toward which they were initially moving. The oscillations are associated with the 
essentially inertial successive recrossings of the higher-energy particles. 

All the arguments developed in the reduced system should have a counterpart in the 
extended scheme. In the energy-diffusion-limited case we have considered small k, large 
uj 2 , and arbitrary 7 [subject to the constraint (0)]. Large ui 2 leads to a narrow harmonic 
potential for y [and consequently y(t) remains small], and small k means weak coupling. 
Since the coupling term ky in Eqs. (^) provides the only energy loss channel for the x 
coordinate, large uj 2 and small k therefore lead to low dissipation. 

Let us now move on to the high dissipation regime. The dashed curve in Fig. |2| mimics 
the exponential memory kernel in the diffusion-limited regime. The kernel T(t) has a high 
initial value and decays essentially monotonically (although we are in the oscillatory regime). 
In the case of the exponential memory kernel (||) this regime results when 7 is large. For 
the exponential memory kernel, the choice of the second parameter, r, further determines 
two different regimes of behavior. If the correlation time r is small, such that 7/r < 1, 
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the transmission coefficient decays monotonically and one is said to be in the non-adiabatic 
regime. This is also the unique high-dissipation behavior associated with the Markovian 
problem. On the other hand, if r is small, such that 7/r > 1, one is in the caging regime 
in which the behavior of the transmission coefficient is oscillatory (but quite differently so 
than in the low dissipation regime). The quasi-exponential dashed memory kernel in Fig. [2] 
corresponds to the non-adiabatic behavior since the ratio of the initial value (about 0.8) to 
the decay time of the kernel (about 4.0) is clearly smaller than unity. In order to have the 
oscillatory memory mimic the non-adiabatic exponential memory case as in the figure we 
require 7 to be small (since now 7 plays the role that 1/r did before), and f2 to be small 
as well (to minimize oscillatory effects). A typical set of values that meets these various 
conditions is u 2 = 0.01, k = 0.75, and 7 = 1.74, which leads to Q = 0.0557. 

The associated transmission coefficient for these parameters exhibits the typical features 
of the non-adiabatic regime, as shown by the monotonic curve in Fig. ||, namely, a smooth 
rather rapid decay to the equilibrium value. As in the case of an exponential memory, this 
decay in the non-adiabatic regime looks Gaussian rather than exponential at short times. 

In the extended system small oj 2 means that the potential in y is very wide and for 
this reason y{t) easily achieves large values. Large k represents strong coupling, and the 
combination of both conditions leads to high dissipation for the x coordinate. 

As noted above, the other regime found for an exponential memory kernel with high 
dissipation is the caging regime, which there occurs when both 7 and r are large, with 7/r > 
1. Interestingly, the constraints on the parameters and the shape of the oscillatory friction 
kernel do not admit this regime. This can be understood from the following argument. 
Caging is achieved when T(t) is essentially constant over some substantial time range so 
that the friction integral in Eq. (|l|) over this time can be approximated as a linear force on 
x(t) and such that the resulting potential becomes monostable. In our case, this resulting 
potential V r (x) would be 

- > 2 - 1)2 + + - l y - 1)2 + ^ < 24 > 

From this expression it is easy to deduce that the potential V r {x) loses its barrier when 
k > 1. The combination of the condition that T(t) behave roughly as a constant for some 
time interval (7 small) and that the resulting potential lose its barrier during this time 
(k > 1) would lead to a caging regime with effective caging potential frequency io cag = 
\/k — 1. However, this combination of conditions can not be satisfied with an oscillatory 
memory. If we increase the value of k above 1, we also have to increase 7 (O has to remain 
small to avoid pronounced oscillations), but this in turn leads to the rapid exponential 
decay of T(t). It is thus not possible to achieve the conditions for the caging regime with 
trigonometric oscillatory friction. The caging regime is easily captured in the hyperbolic 
case (cf. Appendix [A]), since then T(t) can take on a very high initial value that can be 
sustained for a long time. 

We have thus seen that the form of F(t) and the constraints on k, to 2 , and 7 determine 
which regimes typical of exponential memories can also be captured with an oscillatory 
memory. The requirements described so far have been met by either choosing T(0) to be 
small (low dissipation) or large (high dissipation) while minimizing the amplitude of the 
oscillations. 

New behaviors for the dynamics of n(t) may appear in parameter regimes that emphasize 
the oscillatory behavior of T(t). To provide such emphasis we minimize the damping effects 
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of the exponential part by choosing 7 to be small. Further choosing a very small value for 
uj 2 and a medium value for k (large enough to get high amplitudes of oscillation of T(t), but 
limited so that the frequency ft, is not too high) we obtain the oscillatory friction shown as 
the dotted kernel in Fig. ||. Since 7 and to 2 are very small, the frequency of the oscillations 
is 

n ss Vk. (25) 

In this regime at low temperatures an entirely new temporal behavior emerges for K(t). We 
generically call this the "stair-like" regime. It is shown for a typical set of parameter values 
in Fig. H (corresponding to those of the dotted curve in Fig. ^). The main feature of this 
new behavior is that n(t) exhibits a "stair" shape, namely, it decays via a series of steps 
followed by plateaus. The explanation of this behavior, including the period of the steps, 
the dependences on the parameters 7, to 2 , k, and ksT, and the connections with the other 
regimes are presented in detail in Sec. ^. 



5 Approximations 

In this section we lay the groundwork for the arguments invoked in the next section, where 
we discuss the stair-like regime in detail. Our explanations are semi-quantitative, that 
is, we do not develop a theory that reproduces the stair-like curve in the figure in all 
its details. We mention this because in fact such theories are available for the two other 



curves [11, 12, 14, 1{|. In the high dissipation regime the KT theory predicts the monotonic 
decay in the non-adiabatic regime, and this prediction, shown in Fig. |3[ is seen to be 
quantitatively very good (see also jo],]^]). In the low dissipation regime we have shown 
that our theory also leads to very good quantitative agreement with numerical results for 
the entire time evolution of n(t) fljH|T^|. We have not derived such a detailed formula for 
the stair-like regime, but we nevertheless have been able to gain considerable understanding 
of this behavior, and this is what we shall present. 

Our insights turn out to be most complete if we invoke both representations of the 
oscillatory problem, the extended as well as the reduced. Furthermore, an understanding 
of the early time dependence of n(t) in both of these representations turns out to be very 
helpful, even if the approximations that are invoked are not valid over the entire time regime 
- the breakdown of approximations can also yield useful insights. We thus first turn to the 
early time behavior. 

Consider first the reduced representation. KT theory focuses on the way in which 
particles subject to dynamics of the generic form (|l|) with the potential approximated by 
a parabolic barrier diffuse to one side or the other of the barrier. Since their analysis is 
restricted to a parabolic barrier (rather than a bistable potential), the theory is appropriate 
only for high dissipation, that is, when the reaction coordinate never recrosses the barrier 
once it has left the barrier region. In other cases KT theory may (and indeed does) capture 
only the initial decay, typically up to the first plateau value of ft(i), but it does not capture 
the asymptotic values K s t- This is seen in Fig. ||, where the KT theory predictions are 
shown for each of the generic transmission coefficients. KT theory works very well for 
all times for the non-adiabatic high dissipation curve, and captures the initial decay in 
the energy-diffusion-limited (dashed curve) and stair-like (dot-dashed curve) cases. These 
initial agreements are fairly typical for all parameter values. 

Consider now the extended representation (||). We introduce an even simpler approxi- 
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Figure 4: Left panel: The dotted and dashed curves are numerical simulations of the 
transmission coefficient n(t) vs t in the energy-diffusion-limited regime. In both cases u> 2 = 
1.0, k = 0.14 and fcgT = 0.025. Dotted curve: 7 = 0.667 (same value as in Fig. ||); dashed 
curve: 7 = 0.05. The solid line, which reproduces both early time behaviors very closely (see 
inset), is obtained from our early-time approximation for the extended system, Eq. (pq). 
Right panel: the dotted curve is our simulation in the stair-like regime (cf. Fig. ||). The 
solid curve is the early-time approximation. 



mation in this representation that also captures the early time behavior when the dissipa- 
tion of energy in the x-coordinate is slow and that facilitates our analysis of the stair-like 
regime. The approximation is based on three main assumptions, all appropriate only at 
short times. One is akin to the argument we used earlier in the energy-diffusion-limited 
problem, namely, that the main influence of the temperature arises from the initial thermal 
distributions. Thus, as long as the initial distributions are chosen correctly, that is, accord- 
ing to Eqs. (|A6j), ( |A7[ ), (|i~9|), and (|20|), the thermal effects in the form of the explicit random 
force acting on the solvent coordinate can be omitted from the dynamical equations. The 
second is the omission of the dissipation term, i.e., we set 7 to zero. Note that this is the 
dissipative force on the solvent coordinate; the principal initial dissipative channel for the 
reaction coordinate x(t) is its coupling to the y coordinate via k, and this is retained. The 
third is to use a parabolic barrier to approximate the potential. With these assumptions 
the initial decay of the transmission coefficient is due to the low-energy particles (i.e. those 
barely above the barrier) that are pulled by the y coordinate in a direction opposite to the 
one indicated by their initial velocity, as described by the simplified deterministic coupled 
linear equations 

x(t) = (1 — k)x + ky 

y(t) = -(lo 2 + k)y + kx. (26) 
With the initial distributions ( |A6|) , (|A7| ), (|l^) and (p0|) we can then use the form Eq. (^|) 
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for the transmission coefficient to write 



/oo poo poo 

dv yo / dy / dv xo P(v xo ) P{y ) P(vyo)sgn[x(t;v xo ,y ,v yo )} (27) 
-oo J —oo JO 

where sgn[x] is the sign function, that is, sgn[x] = +1 if x > and —1 if a; < 0, and 



x(t;v X o,y ,Vyo) is the solution of Eqs. (26) with initial conditions v xo , y , v yo , and x = 0. 

The left panel in Fig. shows two simulations in the energy-diffusion-limited regime 
along with the results of (pTj) with (|26|). The early time agreement in this regime is clearly 
excellent, as seen in the detail inset. Note that both simulations exhibit the same early-time 
behavior, even though the value of 7 is very different for the two cases (and not particularly 
"small" in one of the two cases). Clearly, in this regime the values of k and oj 2 determine 
the early time behavior of the transmission coefficient. 

More importantly for our purposes here, the right panel of Fig. ^ shows similar early-time 
agreement between the stair-like numerical results and the approximation. The agreement 
extends through the first plateau. Note that the approximation captures the (slightly) 
non-monotonic behavior of the simulation results. The agreement between the two curves 
provides the basis for our analysis of the stair-like regime in the next section. 



6 The Stair-like Regime 

We saw in Sec. |I| that the stair-like regime is achieved when 7 and u> 2 are small and the 
temperature is low. The main feature of this new behavior is that the transmission coefficient 
shows progressive decays connected by plateaus. The length of the plateaus (determined 
by a time period that we call T K ), the depths of the decays, and all the characteristics that 
define this regime depend on the values of the parameters, principally k. Our understanding 
of this regime is based on argumentation that relies mainly on the extended system, although 
some of the arguments can easily be translated to the language of the reduced scheme. 

6.1 Trajectories 

A particularly helpful view of the process is gained by looking at explicit trajectories, 
as illustrated in Fig. ^. The solid trajectory x(t) in the left panel illustrates the typical 
repeated recrossings in the low dissipation energy-diffusion-limited regime. The dotted 
curves correspond to two typical trajectories in the non-adiabatic regime. Trajectories in 
this regime almost never recross the barrier. Those that do recross the barrier do so at short 
times (before straying far from x = 0), and typically do so only once. These trajectories 
reinforce the idea that in this regime particles are quickly trapped in one well or the other 
due to the high dissipation. As seen in the right panel of Fig. ||, the trajectories in the 
stair-like regime are considerably more complex - this complexity distinguishes the stair- 
like dynamics from the other regimes. For example, in this new regime one finds x-particles 
that remain localized over one well (even though they have sufficient energy to cross the 
barrier) and that after circling there several times may suddenly recross the barrier. This 
behavior is not found in any other regime studied so far. 

Studying additional x-trajectories besides those shown explicitly in the right panel of 
Fig. |5] leads to the realization of a number of important points. First, we note that 

• the x-particles cross the barrier only at certain specific times. 
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t t 

Figure 5: Several typical trajectories of x for different regimes in the oscillatory friction 
problem. The parameters are those of Fig. |. Solid curve in left panel: energy-diffusion- 
limited regime, clearly showing repeated recrossings. Dotted curves in left panel: two 
examples of non-adiabatic behavior. Right panel: several trajectories associated with the 
stair-like regime. 

For our typical parameters, these times are t = 10, 16, 25, 31 and 40 (approximately), which 
coincide with the decay times T K for the associated n(t) in Fig. |3|. Moreover, since we are 
working with low temperatures, the energy of the x-particles is typically not large enough 
to recross the barrier many times. Actually, we have observed that 

• most of the x-particles that cross the barrier do it only once. 

Indeed, only about 1% of the particles show multiple recrossings in our typical example. 
The fact that most particles do not return to their original well once they have crossed the 
barrier leads to 

• steps rather than oscillations. 

On the other hand, small friction leads to very slow energy loss, and it is for this reason 
that 

• even at long times crossing the barrier is still possible. 

This combination of features characteristic of the energy-diffusion-limited and non-adiabatic 
regimes ultimately leads to the 

• appearance of successive steps and plateaus. 

At this point there are two obvious questions about this regime: i) Why do we see decays 
only at fairly sharply defined specific times and what are these times? In other words, how 
is the period T K determined? ii) What determines the depth of each decay? The answers to 
these and other questions are given in the following subsections by considering the effects 
of varying the parameters of our model. 

However, before going ahead, we should first understand how an x-particle can be 
trapped in one well in spite of having enough energy to cross the barrier, as well as how an 
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a>particle can cross the barrier seemingly without having enough energy to do it. The rea- 
son for both strange situations is that the behavior of the ^-particles may depend strongly 
on that of their associated y-oscillators. A given x-particle with energy greater than the 
barrier height can be "trapped" in one well because each time it goes toward the barrier its 
coupled oscillator pulls it back. Conversely, a given x-particle with apparently insufficient 
energy can cross the barrier by being pulled by its y oscillator. Thus, the coupling between 
x and y may determine the times at which the particles cross the barrier and therefore the 
times for the decays of n(t). 

6.2 Dependence on k 

To deduce the role of the coupling constant k, we depart from our typical value (k = 0.3) 
and consider the trajectories for two cases on either side of this value but that still essentially 
preserve the stair-like behavior. We call them the large-/c case (k = 0.5) and the small-/c 
case (k = 0.1). The parameter k determines the extent to which x and y particle dynamics 
are in synchrony. 

Since K(t) contains information averaged over ensembles of particles, we are interested 
in the average trajectories of both x and y coordinates. We thus plot (x(t)) and (y(i)), 
where (• • • ) here means an average over all the particles that are in the right well (x > 0) at 
each given time. Fig. |6| shows averaged trajectories for the small-A; (left panel) and large-A; 
(right panel) cases. 




Figure 6: Left panel: {x(t)) (solid curve) and {y(t)) (dotted curve) for the small-fc case. 
Right panel: (x(t)) (solid curve) and (y(t)) (dotted curve) for the large-fc case. 

For the small-A; case we readily observe that the motion of x(t) and y(t) are essentially 
uncorrelated. Each dynamics proceeds with a different principal frequency of oscillation. A 
Fourier analysis of the trajectories reveals that the peak frequency for x(t) is 1.022 while 
that of y(t) is 0.306. These frequencies can be associated with two characteristic frequencies 
of our problem. That of x(t) corresponds to the frequency of the particle in the bistable 
potential, namely 1.022 « 27r/T sem j, where T sem i is roughly the average semiorbit time 
for an ensemble of particles above the barrier in the double- well potential. In Ref. [11| 



we have shown that the semiorbit time for a particle at an energy e above the barrier is 
t e = ln(16/e) + O(elne). An average of this time over a thermal distribution then directly 
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yields T sem i « 3.35 ... — In ksT. The frequency of y(t), on the other hand, coincides with 
the frequency of T(t), namely 0.306 ~ £1 ~ y/k + u; 2 . This frequency characterizes the 
motion of y(i) in the second equation for the extended system, Eq. (0), when the coupling 
contribution is neglected [see also Eq. ([2^)]. The x-particles can cross the barrier with 
frequency 2ir/T sem i since they do not care where their associated y-particles are. Thus in 
this weakly coupled regime 

T K « T semi « 3.35 ... - In A B T. (28) 

The large-A case exhibits a different behavior. In this case x(t) and y(t) are essentially 
synchronized. We see in Fig. |6| that x(t) has two characteristic periods: a shorter one 
associated with motion in the bistable potential (27r/T sem j) and a longer one that matches 
that of y(t). The frequency of the latter is 0.624 and coincides with \fk + uj 2 ^Vk. Thus, 
the motion of x is now dominated by the dynamics of its coupled oscillator. The consequence 
of the strong coupling is that now 




Figure 7: Surface contours for the potential V(x, y) for the small-A; (left panel, k = 0.1) 
and large- A; (right panel, k = 0.5) cases. The coordinates x and y are represented along the 
horizontal and vertical axes respectively. 

These ideas can be further supported by considering the two- variable potential, Eq. (f|), 
drawn in contour form in Fig. [7| for the small-A; (left panel) and large-A (right panel) cases. 
These plots clearly illustrate the correlations between x and y (or the lack thereof). When 
the system is in one of the two two-dimensional wells, x and y remain more tightly bound 
in the large-A case than in the small-A case. In particular, when k is small the y-particle 
can move away from x even when the system has already fallen into one well. 

Further, consider the likely pathways followed by the system as it crosses, say, from the 
right to the left. In the small-A; case the likely path is for y to decrease first (perhaps even 
to negative values), followed by a change of x from x > to x < 0. On the other hand, 
in the large-A; case it is easier for x to first move from x > to x < to be followed by 
y. Therefore, in the large-A; case the crossing rate is determined by the frequency of y(t) so 
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that T K ~ 2ir/yk; in the small-A; case the crossing rate is limited by the motion of x(t) and 
hence T K ~ T sem i. 

Figure || shows the time-dependent transmission coefficient for the two cases. The times 
T K 7.55 (small-A; case) and T K m 10 (large-A; case) obtained from the above arguments 
are consistent with the steps in the figure (measured from mid-point to mid-point), most 
clearly in the length of the first step. 




0.0 10.0 20.0 30.0 40.0 50.0 0.0 10.0 20.0 30.0 40.0 50.0 

t t 



Figure 8: Left panel: n(t) for the small-A; and large-A; cases discussed in the text. Right 
panel: Transmission coefficient for the small-A; case for three different barrier heights. The 
solid curve is the same as the small-A; case in the left panel. 

In addition to the step period differences, the left panel in Fig. ||| illustrates the k- 
dependence of the depths of the steps in the stair-like transmission coefficient. The steps 
are clearly deeper when k is small. The reason for this is that the total system loses energy 
by dissipation only through the y coordinate. Although both cases in the figure correspond 
to the same value of 7, the x-particles can retain energy for a longer time when k is small. 
This allows more x-particles to cross the barrier, and it allows them to do so at later times. 
The deeper and more numerous clear steps in the small-A; case are a direct manifestation of 
these features. 

Two further points should be noted. One is the symmetry of the semiorbit time t e with 
respect to s. That is, the semiorbit time of a particle with an energy e above the barrier is 
the same as the orbit time of a particle with energy e below the barrier (for small s). This 
symmetry is important because it allows particles to remain in synchrony; otherwise the 
steps in n(t) would be blurred. The other point is the dependence of n(t) and consequently 

— 1/2 

of the period T K on barrier height. In general, t e ~ V ln(16Vo/e) (which reduces to 
our previous expression when Vq = 1), and an average of t e over a thermal distribution of 
particles above the barrier yields the generalization of Eq. (p§|) 

' (In I6V0 + 0.5772... -In k B T) . (30) 



^0 



The right panel in Fig. || shows the transmission coefficient in the small-A; case for three 
values of the barrier height. The corresponding period estimates for T k (Vq) obtained from 
Eq. pi) are T K (0.5) = 9.70, T«(1.0) = 7.55, and T K (2.0) = 5.83. These decreasing periods 
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with increasing barrier height are clearly consistent with the numerical results. 
6.3 Dependence on 7 

We have seen that successive steps in the transmission coefficient arise because the particles 
lose their energy slowly. This requires 7 to be small - but not too small (cf. below). Indeed, 
if 7 is decreased we expect particles to lose their energy even more slowly, which leads to a 
larger number of deeper steps. However, as 7 continues to decrease we expect to begin to 
see particles that cross the barrier more than once before becoming trapped. This leads to 
oscillations in n(t) and, eventually, to energy-diffusion-limited behavior. Deeper steps and 
the first appearance of small oscillations with decreasing 7 are clearly evident in the left 
panel of Fig.[| 

Conversely, if 7 is increased, particles lose their energy more rapidly, and fewer particles 
cross the barrier at all; those that do so cross at most once. In this case, as seen in Fig. pi 
the steps are less deep and almost disappear at long times. Indeed, the limit of the stair- 
like behavior with increasing 7 is the monotonic non-adiabatic regime, where only a few 
particles cross the barrier and they do so at very short times. 




Figure 9: Left panel: dependence of the stair-like transmission coefficient on the dissipation 
parameter 7. The parameters are otherwise those of our typical example in Fig. |3|: cj 2 = 
0.01, k = 0.3, ksT = 0.015. Dotted curve: 7 = 0.05; solid curve: 7 = 0.01; dashed curve: 
7 = 0.1. Right panel: dependence of the stair-like transmission coefficient on temperature; 
the other parameters are those of Fig. ||. 



6.4 Dependence on ksT 

Finally, we consider the temperature dependence of the transmission coefficient in the stair- 
like regime. As temperature is increased, all else remaining the same, there is a greater 
number of more energetic a>particles above the barrier. Two aspects of their behavior 
dominate the resulting transmission coefficient. One is that the particles now have a greater 
range of semiorbit times t e ; the other, more important, effect is that particles are now 
sufficiently energetic that they can recross the barrier more than once. These are precisely 
the features that lead to the typical oscillatory behavior of the transmission coefficient 
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in the energy-diffusion-limited regime, and it is towards this behavior that the stair-like 
regime tends with increasing temperature. The right panel in Fig. |9] shows this progression 
very clearly: the highest temperature results look very much like the earlier curves for the 
energy-diffusion-limited case. We should note the deeper first decay in the /cbT = 0.5 curve 
than observed in our earlier illustrations. This is due to the fact that here we have chosen 
uj 2 to be very small (a requirement for the stair-like regime). This causes the initial thermal 
distribution of y(0) to pull back x-particles more effectively than in our earlier example, 
and this in turn leads to the deeper decay. 

6.5 Arguments in the Reduced System 

Since the extended [Eq. (g)] and reduced [Eq. (|])] systems are entirely equivalent, it is of 
course possible to explain the new stair-like regime in terms of quantities and equations 
associated with the reduced representation. It is perhaps somewhat more cumbersome 
and less transparent, but the extended representation analysis offers a helpful guide. For 
example, it is useful to realize that a "negative friction" [i.e., a negative value of the memory 
kernel in Eq. (|l|)] in the reduced representation is associated with the situation where in 
the extended system the y-coordinate pulls the x-particle in the direction of x. 

To explain the different decay periods T K in the reduces representation we note that the 
memory kernel T(t) is proportional to k 2 and the random force F(t) is proportional to k. 
To cross the barrier, an x-particle must be moving towards x = 0. When k is small, the 
dynamics of x as it moves in the barrier region is thus dominated by the bistable potential 
V e ff(x). In particular, the decay periods in K,(t) are determined by the frequency of the 
particles moving in the bistable potential with energies slightly larger and slightly smaller 
than the barrier height, which directly leads again to the earlier estimate T K = T sem i [cf. 
Eq. (|28|)1. The slow dissipation of x-energy associated with small k allows for many deep 
steps in K,(t). 

As k increases, the bistable potential becomes relatively less important and the first and 
third terms on the right of Eq. ([l]) increasingly dominate the dynamics of x(t). The steps 
in n(t) then acquire the period T K ~ 2ir/n ~ 2ir/ 1 \pk associated with the friction kernel. 
The more rapid dissipation of x-energy associated with larger k leads to a small number of 
shallow steps. 

The dependence on 7 in this representations is quite clear. When 7 increases, the mem- 
ory kernel decays more rapidly and the oscillations in T(t) become irrelevant, thus leading 
to non-adiabatic behavior of n{t). Decreasing 7, on the other hand, leads to pronounced 
oscillations (and at times negative values) of T(t). As a result, even particles that start out 
with energies too low to cross the barrier early may do so at a later time, thus explaining 
the step structure of K(t). The temperature dependence can also be understood: for a 
given (low) 7, increasing the temperature leads to a greater number of particles above the 
barrier that can recross more than once before becoming trapped. The steps then become 
oscillations and the energy-diffusion-limited behavior is recovered. 

7 Conclusions 

In this work we have analyzed the time dependent transmission coefficient for the capture 
of a particle in one or the other well of a bistable potential as described by the generalized 
Kramers equation Eq. (|]) with an oscillatory memory kernel. The time dependence of the 
transmission coefficient depends sensitively on the parameters of the model. The equivalence 
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of this model to an "extended system" wherein the reaction coordinate is linearly coupled 
to a nonreactive coordinate which is in turn coupled to a heat bath, Eq. (||), facilitates the 
understanding of the various time dependences that are observed. 

The different behaviors observed for the transmission coefficient in various parame- 
ter regimes are summarized in Fig. |3| The non-adiabatic (monotonic decay) and energy- 
diffusion-limited (oscillatory decay) behaviors have been encountered earlier in the classic 
Kramers problem |IT| , |il^ , |i4| ,|l6| and in the generalized Kramers problem with an exponential 
memory kernel [12-^4, 16|. The non-adiabatic decay is observed when the reaction coordi- 
nate loses its energy rapidly, so that particles cross the barrier only at early times and at 
most once before becoming trapped. The oscillatory behavior is observed when the reac- 
tion coordinate loses its energy slowly, thus allowing several recrossings of the barrier before 
trapping. These regimes are associated with parameter values that suppress the oscillations 
of the memory kernel. There is a third behavior observed with exponential friction, the 
caging regime, which is not observed with an oscillatory memory friction. 

The third behavior shown in Fig. ^, which consists of a stair-like decay of the transmis- 
sion coefficient, is peculiar to the oscillatory memory friction and occurs when the oscilla- 
tions in the memory kernel are pronounced. This behavior is observed when particles cross 
the barrier at most once, but not necessarily at early times. In turn, this can be explained 
by the fact that particles that at one time may not have enough energy to cross the barrier 
may acquire sufficient energy to do so later via their coupling to the nonreactive coordinate 
(or, equivalently, when the oscillations in the memory kernel periodically lead to negative 
values of the kernel). Although the particles cross the barrier at most once, and not nec- 
essarily at early times, the crossing events can only occur at fairly sharply defined time 
intervals that we call T R . Hence the appearance of fairly sharp steps in the transmission 
coefficient. We explain in detail the conditions that lead to the stair-like behavior, the way 
in which the step time T K and the step depths depend on the parameters of the system, and 
the way in which this behavior tends to the energy-diffusion-limited or non-adiabatic cases 
as parameters are modified. 

If there were no barrier crossings at all in the Kramers problem, the transmission coef- 
ficient would be unity. Single barrier crossings only at early times lead to monotonic decay 
of the transmission coefficient. Single recrossings that are possible only at specified time 
intervals T R lead to the new stair-like regime. Multiple recrossings lead to oscillatory be- 
havior. The time dependence of the transmission coefficient clearly provides an interesting 
mirror for the barrier crossing dynamics of the generalized Kramers problem. 
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A Extended vs Reduced Model 

Here we present the analytical details connecting Eqs. (|l]) and (^). Formal solution of 
Eq. (|) gives 

= Vyo ~ yoA 2 gAlt + VoXl ~ V yo e x 2t 

Ai — Xo Ai — Ao 



+ 



Ai - A 



2 7o 



di' ( e M*-*') - e *»(*-0) + /(*')], (AI) 



where the first two terms on the right hand side correspond to the homogeneous solution that 
depends on the initial conditions y(0) = y a and y(0) = v yo . The last two terms correspond 
to the inhomogeneous solution. The term proportional to k leads to the memory friction 
term and the term containing f(t') is associated with the colored noise in the reduced model. 
The roots A,- are 



Vl.2 



7 , //7 n2 



± ± _J2_ km (A2) 



Substitution of this formal solution in Eq. (^) and regrouping of terms directly leads to 
the reduced model (|l|) with the memory kernel 

^,2 / e Ait e A 2*\ 
A2 — Ai \ Ai A2 / 

We also get the explicit form for the effective potential of the reaction coordinate, 

V eff {x)=V{ X ) + \-^x\ (A4) 

However, the following extra initial conditions must be fulfilled in order to avoid transient 
terms in the reduced model: 

kuT 

x(0) = ; < v yo >=< y >= ; < v 2 >= k B T ; < y 2 >= (A5) 

The brackets here indicate averages over initial distributions. The following initial distri- 
butions for the solvent coordinate are consistent with these requirements: 



and 



1 / vl 



P(v vo ) = , = exp I . (A7) 

1 y ' y«^T ^ \ 2k B T J V > 

The explicit reduction (integration) thus readily leads to the observation that the initial 
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conditions for y are the thermalized solutions of the homogeneous differential equation 

y + jy + (w 2 + k)y = 0. (A8) 

That is, we must thermalize the solvent coordinate evolving in the combined intrinsic and 
coupling potential with the reaction coordinate fixed at x = 0. 

At this point, a distinction should be made between the following two behaviors of the 
friction kernel. The first is the underdamped case, where the condition 

n 2 = to 2 + k- (J) 2 > (A9) 

leads to complex values for Ai and A2, 

Ai )2 = -|±ifi, (A10) 

which in turn leads to a trigonometric form for the memory kernel (sometimes called the 
"trigonometric case" ) : 

T(t) = fc e'l* sin Qt + cos fit) . (All) 
The second behavior, the overdamped case, results when 

A 2 = -u? - k > 0. (A12) 

In this case the values of Ai and A2 are real, 

Ai, 2 = -|±A, (A13) 

and therefore the memory kernel has a hyperbolic form (sometimes called the "hyperbolic 

case"): 

I7t) = k - e~2* f^- sinh At + cosh At) . (A14) 
uj 2 + k V2A / 
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